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ABSTRACT 

The  recent  and  imminent  launch  of  the  SMOS  and  Aquarius 
satellites  carrying  microwave  L-band  radiometers  provides 
an  opportunity  to  map  Sea  Surface  Salinity  (SSS)  globally 
with  an  expected  error  <  0.2  psu.  However,  the  accuracy  of 
retrieved  SSS  depends  critically  on  brightness  temperature 
(Tb)  corrections  for  sea  surface  roughness  (SSR)  effects. 
This  paper  assesses  the  performance  of  representative 
roughness  correction  models  when  eompiircd  with 
published  data,  and  applied  to  recently-acquired  airborne  L- 
band  radiometer  data.  One  type  of  model  currently  being 
used  to  process  SMOS  data  combines  a  wind-driven  gravity 
wave  spectrum  that  describes  SSR,  with  an  electromagnetic 
(EM)  model  that  determines  microwave  emissivity,  to 
predict  the  Tb  roughness  increment  relative  to  the  flat  sea 
response.  We  find  that  selection  of  both  the  spectral  and 
emissivity  models  strongly  influences  the  resulting  (~1  K) 
Tb  errors.  We  conclude  that  more  accurate  modeling  of 
short  wavelength  spectral  components  and  their  EM 
influence  is  needed,  to  reduce  these  errors  to  acceptable 
levels. 

Index  Terms —  L-band  Radiometer,  Sea  Surface 
Salinity,  Emissivity,  Surface  Roughness,  Wave  Spectra 

I.  INTRODUCTION 

The  launches  of  ESA's  Soil  Moisture  and  Ocean  Salinity, 
SMOS,  L-band  microwave  satellite  [8]  in  Nov.,  2009,  and 
NASA's  Aquarius  satellite  [19]  in  late  2010,  promise 
monthly  global  maps  of  Sea  Surface  Salinity  (SSS)  over  the 
deep  ocean  with  0.2  psu  precision  at  100  km  resolution. 
However,  SSS  retrieval  accuracy  depends  critically  on 
corrections  for  SSR  effects.  These  effects  change  the  L- 
band  (~21  cm  wavelength)  sea  surface  brightness 
temperature,  Tb,  predicted  by  'flat  sea'  emissivity  models 
[17],  as  much  as  do  open  ocean  surface  salinity  variations. 
Thus,  they  require  careful  correction.  This  paper  describes 
the  evaluation  of  roughness  correction  models  designed  for 


retrieving  SSS  from  L-band  microwave  radiometer  Tb 
measurements  over  rough  seas.  The  paper  first  describes  the 
airborne  instrumentation,  and  the  wave  spectrum  and 
emissivity  models  used.  Then  results  from  several  models 
arc  presented  and  compared  using  published  data  and  recent 
NRL  Salinity,  Temperature  and  Roughness  Remote 
Scanner,  STARRS,  airborne  salinity  mapper  campaign 
measurements.  It  is  shown  that  use  of  Hwang's  recently 
published  wave  spectrum  [12]  to  drive  the  SSA/SPM  [15] 
[16]  [22]  [25]  emissivity  model  can  eliminate  the  need  for 
an  overall  factor  of  two  increase  in  spectral  intensity  used  in 
previous  roughness  correction  models. 

2,  METHODS 

2.1.  Airborne  Instrumentation 

STARRS  comprises  L  and  C-band  microwave  radiometers, 
an  infrared  (IR)  radiometer,  integrated  GPS  receiver,  and 
gyro.  The  combined  measurements  are  used  to  retrieve  sea 
surface  salinity  (SSS),  temperature  (SST)  and  roughness. 
The  L-band  radiometer,  used  primarily  to  retrieve  SSS,  is  a 
multi-beam  system  sensing  microwave  thermal  emission 
from  the  sea  surface  within  a  24  MHz  wide  protected  band, 
centered  at  a  frequency  of  1.4  GHz.  The  6  antenna  bcam.s, 
with  15  deg  half  power  beam  width,  point  dow  nward  and  to 
cither  side  of  the  aircraft  at  incidence  angles  of +/-  7.  22  and 
38  degrees.  The  nadir- viewing  IR  radiometer  senses  SST 
from  thermal  emission  in  the  8-14  and  9.6-11.5  mieron 
bands.  For  typical  aircraft  speeds  of  80  m/s  and  altitudes  of 
2600  m,  the  beam  geometry  yields  a  5.2  km  swath  width. 
Sec  [I]  and  [20]  for  other  instrument  and  processing  details. 

2.2.  Roughness  Emissivity  Models 

Various  emissivity  models  are  presently  available  to  correct 
for  the  adverse  clTccts  of  roughness-enhanced  emission  on 
microwave  SSS  retrievals.  These  include  rigorous  [21], 
asymptotic  [16]  [22]  [26]  and  empirical  model  types  [2]  [3] 
[9].  Three  models  of  the  last  two  types  arc  implemented  in 
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no  (nil) 

Kig  1.  SPM/SSA  calculations  from  dilTcrent  wind-wave 
Spectra.  (Top)  Curvature  spectra  B(k)  of  Elfouhaily  el  al.  [5] 
and  Kudryavtsev  et  al.  [13]  versus  wavenumber  k(rad^m). 
(Bottom)  Predicted  Th  (K)  versus  wind  speed  t/(m/s)  at  10  m 
from  SSA/SPM  for  Inc.Ang.  37  deg.,  Ts  298  K,  S  35  psu, 
Inv.wave  age=0.84. 

the  SMOS  Level  2  processor.  The  asymptotic  models  are 
driven  by  auxiliary  wind  data  and  based  on  a  specified 
wind-wave  spectrum.  Their  accuracy  is  strongly  influenced 
by  the  choice  of  spectrum.  The  empirical  models  arc  driven 
directly  by  auxiliary  wind  or  sea  state  data,  and  calibrated 
using  in  situ  observations.  They  take  advantage  of  SMOS’s 
multi-angle  view  capability  [8].  In  contrast,  Aquarius  will 
model  roughness  empirically  using  radar  cross  sections 
observed  by  an  on-board  L-band  scatteromeicr  [19].  The 
variety  of  models  adopted,  even  for  the  single  mission, 
SMOS,  is  a  reflection  of  recent  issues  and  innovations  in 
this  field,  and  of  the  dominant  role  of  sea  surface  roughness 
corrections  in  the  L-band  radiometer  error  budgets. 

2.3.  Wave  Spectral  Models 

The  asymptotic  cmissiviiy  models  are  combined  with  a 
wind-driven  SSR  spectrum  to  form  a  complete  roughness 
correction  model.  From  the  perspective  of  determining 
microwave  cmissivity,  the  gravity  wave  spectra  remain 
poorly  defined.  The  difficulties  of  specifying  an  optimal 
SSR  spectrum  is  highlighted  in  the  comparison  between 
field  measurements  and  analytical  computations.  For 
example,  Yuch  [26]  selects  the  Durden  and  Vesecky 
spectrum  [6]  for  his  two-scale  model,  but  must  double  the 
roughness  spectral  densities  to  obtain  reasonable  agreement 
with  field  measurements.  Camps  et  al.  [2004]  used  the  [6] 
and  Elfouhaily  et  al.  [7]  spectral  models,  also  doubled,  to 
make  the  calculated  rate  of  change  of  brightness 
temperature  with  wind  speed  agree  roughly  with  their 


Fig  2.  Average  rate  of  change  of  Tb  with  respect  to  wind 
speed,  97b/3UJ0  (K/m/s)  versus  incid.  angle  at  1.41  GHz 
calculated  using  four  different  wind-wave  spectra  for  wind 
speed  range  2  to  14  m/s.  Hollinger  [10]  and  Camps  et  al.  [3] 
field  data  are  superimposed. 

measurements.  While  [6]  produces  mean  square  slope 
(MSS)  similar  to  the  classical  sun  glitter  measurements  of 
Cox  and  Munk  [4],  Yuch  [26]  argued  that  other  well  known 
spectral  models  produce  integrated  total  MSS  about  twice 
the  magnitude  of  the  CM  data,  which  could  have 
underestimated  surface  slopes  that  were  very  steep,  but 
infrequent  [23].  Similar  difficulty  in  specifying  the  SSR 
spectrum  is  also  encountered  in  active  radar  scattering 
computations  [13]  [14]. 

3.  RESULTS 

3.1.  Roughness  Correction  Model  Assessment 

To  investigate  the  importance  of  model  choice  in  computing 
roughness  corrections  for  SSS  retrieval,  we  compared 
results  from  the  Two-Scalc  Model  (TSM)  of  ^'ueh  [26]  and 
RcuFs  (pers.  comm.)  code  for  the  SPM/SSA  asymptotic 
model  [16][22],  both  of  which  arc  implemented  in  the 
SMOS  L2  processor.  These  models  also  employ  different 
wind-wave  spectra  to  describe  sea  roughness.  SPM/SSA  is 
driven  by  the  Kudryavtsev,  cl  al.  spectrum  [  1 8]  and  TSM  by 
the  Durdcn-Vcsceky  spectrum  [6]  (multiplied  by  two,  as 
discussed  above),  with  Gaussian-distributed  long  wave 
slopes  [26].  A  significant  difference  in  the  Tb’s  predicted  by 
TSM  versus  the  SPM/SSA  model  appeared,  particularly  for 
H-Pol.  At  50  deg.  incidence  angle,  TSM  predicts  -  2  K 
lower  Tb  influence  than  SPM/SSA  (-4  psu  SSS  error!). 
Such  under-prediction  by  TSM  has  previously  been 
reported  by  other  investigators. 
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Fig  3.  (Top)  STARRS  SST  transect  from  Chesapeake  Bay 
entrance  to  mid  Gulf  Stream.  (Bottom)  MODIS  satellite  SST 
image  showing  transect  crossing  Gulf  Stream  (Red  line  is 
STARRS  transact,  black  line  is  a  satellite  data  dropout). 

3.2.  Effect  of  Spectral  Model 

To  illustrate  the  effect  of  input  spectrum  choice  for  a  single 
asymptotic  emissivity  model,  we  show  two  surface 
curvature  spectra  (Fig.  1)  and  resulting  predictions  of  L- 
band  H-Pol  Tb’s,  derived  using  ReuFs  implementation  of 
the  SPM/SSA  model  [22].  The  input  wave  spectra  were 
those  of  Kudryavtsev  et  al.  (K)  [18]  and  Elfouhaily  ct  al.  (E) 
[7],  The  resulting  Tb  predictions  differ  by  -1  K  for  an  SSS 
error  of  about  2  psu  under  typical  temperate  conditions  (V- 
Pol  errors,  not  shown,  were  of  similar  magnitude). 
Comparison  of  Tb  wind  sensitivity,  which  results  from 
using  E  and  K  along  with  the  Donclan  (D)  [5]  and  Hwang 
(H)  [12]  spectrum  to  drive  SPM/SSA  (Fig.  2),  shows  that  H 
perfonns  best  for  H-Pol  Tb  and  is  competitive  with  K  for  V- 
Pol.  considering  field  data  spread.  Thus  H  provides  the  best 
overall  perfonnance,  while  both  models  produce  wind 
sensitivities  appreciably  larger  than  those  predicted  using 
the  E  spectrum. 

3.3.  Application  to  Airborne  Measurements 

A  STARRS  transect  crossing  the  continental  shelf  and 
western  half  of  the  Gulf  Stream  conducted  during  NRL’s 
VIRGO  experimental  campaign  in  Dee.,  2006  shows  SST’s 
obtained  from  the  STARRS  and  MODIS  IR  radiometers  and 
Gulf  Stream  location  (Fig.  3),  and  the  corresponding 


Virgo  12  Dec  06  Estimated  Roughness 
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Fig  4.  (Top)  STARRS  SSS  transect  passing  NDBC  buoys 
(red  dots,  W  to  E  CBBV2,  CHLV2  and  440 1 4).  (Bottom)  Bar 
chart  showing  L-band  Tb  [K]  corrections  from  TSM  (blue), 

WS  (green),  and  SPM/SSA  (red)  for  the  buoy  locations,  and 
observed  wind  speeds  [m/s]  (sec  labels). 

STARRS  SSS  transect  (Fig.  4).  The  bar  chart  shows  the  V- 
Pol  Tb  corrections  (Delta  Tb)  computed  from  I  SM  [26],  the 
empirical  WISE  emissivity  model  [2]  (WS),  and  SPM/SSA 
(with  input  K  spectrum),  at  locations  near  NOAA  data 
buoys.  The  wind  speeds  observed  by  the  buoys  during  the 
night  arc  also  shown.  The  corresponding  SSS  correction  in 
this  temperate  region,  in  psu,  is  approximately  2  x  Delta  Tb 
(K).  TSM  consistently  under-prediets  Delta  Tb  with  respect 
to  SPM/SSA  across  the  whole  range  by  a  factor  of  about  2. 
Implied  model-dependent  errors  are  of  order  0.75  K  in  this 
case,  corresponding  to  an  SSS  correction  error  of  about  1. 5 
psu.  This  error  can  be  compared  with  an  SSS  difference  of 
-5  psu  observed  across  the  Gulf  Stream  in  1 099,  using  the 
PALS  L-band  radiometer  [24]. 


4.  SUMMARY 


The  dominant  roughness  inlluence  on  SSS  retrieval  comes 
from  the  shorter  Bragg-scale  components  of  the  wind-wave 
spectrum  (Icm-lm  al  L-band).  Since  these  waves,  in 
particular,  arc  poorly  represented  in  traditional  spectrum 
models,  recent  spectral  model  improvements  have  focused 
on  better  accounting  for  their  effects  using  more 
sophisticated  physical  and  empirical  modeling  [1 2]  [1 8]. 
This  seems  preferable  to  simply  doubling  the  spectral 
intensities  to  predict  roughness  emissivity  inlluence.  As  the 
short  wave  components  arc  modulated  by  the  long  waves. 
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they  are  also  influenced  by  swell  [11],  which  future 
roughness  correction  models  should  also  take  into  account. 

We  conclude  that  it  is  vital  to  choose  the  roughness 
correction  model  and  input  forcing  function  (wind  speed, 
and/or  spectrum)  carefully  to  minimize  this  major  error 
source  for  L-band  SSS  retrieval.  Further  enhancements  to 
the  H  spectrum  wind  dependency  are  being  investigated, 
and  a  rigorous  FDTD  reference  model  is  currently  under 
development  as  an  aid  in  testing  and  enhancing  candidate 
operational  emissivity  models.  The  results  will  be  shared 
with  the  ESA  SMOS  and  NASA  Aquarius  science  teams. 
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